A derivative formula for the free energy function 
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Abstract 

We consider bond percolation on the TI^ lattice. Let M„ be the number of open clusters in 
B(ji) = [— n,n]'^. It is well known that EpMn/ {2n + 1)'^ converges to the free energy function 
k{p) at the zero field. In this paper, we show that (Tp(M„)/(2n +1)'' converges to — (p^(l —p) + 

1 Introduction and statement of results. 

Consider bond percolation on the Z'^ lattice, in which bonds are independently open with probability 
p and closed with probability 1 — p. The corresponding probability measure on the configurations 
of open and closed bonds is denoted by Pp. We also denote by Ep{X) and cTp(X) the expectation 
and the variance of X with respect to Pp. The open cluster of the vertex x, C(x), consists of all 
vertices that are connected to x by an open path. Here an open path from u to u is a sequence 
{vQ,hQ,vi, ...,Vi,hi,Vi+i, ...,Vn) with distinct vertices Vi {Q < i < n) and open bonds hi adjacent Vi 
and Vj+i such that vq = u and = v. For vertex set A, denotes the bonds with both vertices in 
A. Also, 1^1 denotes the cardinality of A, and |Ae| denotes the number of bonds in A^,. We choose 
as the origin. The percolation probability is 

e{p) = Pp{\cm = oo), 

and the critical probability is 

Pc = sup{p : 9{p) = 0}. 
We denote the open cluster distribution by 

enip) = Pp{\C{0)\=n). 

By analogy with the Ising model, we introduce the magnetization function as 

oo 

M{p, h) = l-Y, ^n(p)e""'' for h>0. 

n=0 

By setting /i = in the magnetization function, 

M(p,0) = 9{p). 
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Using term by term differentiation, we also have 

J_im = E,{\Cm; |C(0)| < oo) = x^{p). 

ip) is called the mean cluster size. The free energy F{p, h) is defined by 

oo 1 

F{p, h) = h{l - Ooip)) + J2 -^n{p)e-^'' for /i > 0. 



n=l 



If we differentiate with respect to /i, then we find 

dF{p,h) 



M{p,h). 



dh 

For h > 0, the free energy is infinitely differentiable with respect to p. If /i = 0, F{p, 0) is called 
the zero-field free energy. The zero-field free energy F{p, 0) is a more interesting and more difficult 
object of study since it is believed that there is a singularity point at Pc- By our definition, 

F(p,0) = i?(|C(O)r^|C(O)| >0). 

Grimmett (1981) discovered that the zero-field free energy also coincides with the number of open 
clusters per vertex. Let us define the number of open clusters per vertex as follows. We denote 
by Mn the number of open clusters in B{n). By a standard ergodic theorem (see Theorem 4.2 in 
Grimmett (1989)), the limit 

lim , „ , , , Mn = lim r:;Mn = k(p) a.s. and Li (1-1) 

n^oo\B{n)\ n-^oo (2n + 1)'^ n \yj \ J 

exists for all < p < 1. K{p) is called the number of open clusters per vertex. Grimmett (1981) 
proved that 

Kip) = F{p,0). (1.2) 

k{p), as a function of p, is analytic for p ^ Pc and differentiable on [0, 1] (see Kesten (1982)). 
In particular, k{p) is proved (see Kesten (1982)) to be twice differentiable at Pc for the square 
lattice. In general, physicists believe that the zero-field free energy is twice, but not three times 
differentiable at Pc. 

On the other hand, for k{p), as a limit of random variables, Zhang (2001) showed the following 
central limit theorem. For p G [0, 1], 



Mn - EpMn 



a standard normal distribution. (1.3) 



CTpiMn) 

Zhang also showed large deviations for M„. In this paper, we show another property for k{p). 

Theorem. For < p < 1, 

lim al{Mn)/(.2n + lf = -{pHl-p) +p{l-pfy{p). 
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Remark. If bo is the bond with vertices and (1,0, • • • ,0), and G{b) is the event that there 
does not exist an open path in Zg \ 6 connecting the two vertices of bond b, then by (2.9) in the 
proof of the theorem, we have 

k'(p) = -dPp{g{bo)). (1.5) 
By using Sykes and Essam's formula (1964), we know that for the square lattice, 

k'(0.5) = -1. 

Thus 

lim al^{Mn)/{2n + l)'^ = 0.25. (1.6) 



2 Proof of theorem. 

For any bond 6, let vi{b) and V2{b) be the two vertices of b. Given p, we start by taking the 
derivative of Ep{Mn). Note that 

oo 

EpMn = Y,PpiMn>l), 

1=1 

and the event {Mn > 1} is decreasing. Let {Af„ > l}{b) be the event that 6 is a pivotal bond for 
{Mji > I}- By Russo's formula, note that Pp{Mn > 1) = 1, so 



dEpMn 



dp 



Let 



= -E E Pvmn>m) 

1=2 feeBe(n) 

oo 

= — ^ ^ Pp{Mn = Z — 1 or Hf 6 is open or closed) 

b&Be{n) 1=2 

/oo \ 

= — ^ -^p U = / — 1 or Hf 6 is open or closed . 

b&Be{n) \l=2 J 



£n{b) = {6 is a pivotal bond for the open connection of vi{b) and V2{b) in B{n)}. 



(2.1) 



In other words, if b is open, then fi(6) and V2{b) are connected by open paths. Conversely, if b is 
closed, then vi{b) and ^2(6) are not connected by open paths. Thus, for each b G Be{n), 

Pp M„ = Z - 1 or Hf 6 is open or closed^ = Pp(f„(6)). (2.2) 

Let Qnib) be the event that there does not exist an open path connecting ^1(6) to V2{b) inside 
Be{n) \ b. Then we would have 

Sn{b) = Qn{h). (2.3) 
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-n, n 



Figure 1: The figure shows that Qn{b) occurs, but Q{b) does not occur. The two vertices of b are 

connected by open paths, but they have to reach to the boundary of B{n) before connecting with 
each other. The dotted path indicates the closed bonds that block the connection of open paths 
from vi{b) and ^2(6) inside B(n). 



To see (2.3), if there were such a path, then v\{b) and V2{b) would always be connected by an open 
path whenever b is open or closed. So b would not be a pivotal bond. On the other hand, if there 
docs not exist an open path connecting vi{b) and V2{b) in Bf>{n) \ b, then b should be a pivotal 
bond for the open connection of vi{b) and V2{b). With these observations, 



dEpMn 
dp 



E Ppmb)) = - E PpiOib))- E Pp{On{b)ng^{b)). 



By translation invariance, the ratio of the first term above is 

E Pp{gm/\Be{n)\ = Ppigibo)). 



(2.4) 



We use the following lemma to estimate the second term. 
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Lemma. 

lim V Pp (Qnib) n g^ib)) = uniformly on [0, 1]. 

n^oo \Be[n)\ , , ^ ' 

Proof. Let A{n, m) = B{n)\B{m) for m < n, and let A{n, m) be the closure of A{n^ m). Then 



4^ E p,(a„,.)nac 



< Ib;^ E i'p(5.('-)n5«((,)) + i^ E i'p(5nWn5«(6)) 

< E ^'p(anWna^(6))+0(n-"-^), (2.5) 

where we may assume that n — ^fn is an integer; otherwise we may use \n — ^/n ] to replace n — y/n. 
Let us estimate the first term in the above inequality. For 6 G Be{n — ^/n), £n{b) fl £^ {h) implies 
that there exists an open path from vi{b) and V2{b) without using b, but open paths cannot stay 
inside B{n). In other words, the open path adjacent to vi{b) has to reach the boundary of B{n) 
before reaching V2{b) (see Fig. 1). Similarly, the open path adjacent to V2{b) has to reach the 
boundary of B{n) before reaching vi{b) (see Fig. 1). Therefore, there exist two disjoint open paths 
from Vi{b) and V2{b) such that both reach to the boundary of vi{b) + B{'^ — 1). Let P(6, ^Jn) be 
the event. By using the estimate of Theorem 6.1 in Grimmett (1989), for all p G [0, 1], there exist 
constant C = C{d) and 5 = 6{d) < 0.5 such that 

Pp {V{b, < Cn-\ (2.6) 

By (2.5) and (2.6), 

E Pp{G{b)ng^{b))<0{n-'). (2.7) 



Therefore, the lemma follows from (2.7). □ 

It follows from (2.3), the lemma, and (2.4) that 

^M^ = -{\B,in)\/\Bin)\) E Pp{G{b))/\B,{n)\ + 0{n-'). (2.8) 

Note that 

\Be{n)\ = 2dn{2n + if''^ and \B{n)\ = (2n + l)'', 
so if we let n — >■ oo in (2.8), 

/^'(p) = -Jim E Pp{Gnib))/{2n + lf = -dPpiG{bo)). (2.9) 

b&Bein) 
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Now we estimate the variance of M„. We list the bonds in Be{n) in some order: 

{bi,---,h}- 

We define the independent Bernouni-random variables {oj^bi)} for 1 < i < A; to be uj{bi) = (open) 
or uj{bi) = 1 (closed) with probability p oi 1 — p. Now we construct the following filtration: 

To = {0, CI} C Ti = {cr— field generated by oj{ei)} C, - ■ ■ ,C Tk = {cr— field generated by uj{ei), ■ ■ ■ ,uj{ek)}- 

The martingale representation of M„ — EpMn is 

fe-i 

i=0 

Let the martingale difference be 

Ai^k = [Ep{Mn\Ti+i) - Ep{Mn\Ti)]. 

The variance is 

k-l 
i=0 

Both Mn and Aj^^ can be viewed as functions on [0, 1]^' and [0, 1]*, respectively. So we can write 
-^n(ci, ■ ■ ■ , Cfe) and Aj^fc(ci, • • ■ , Cj) for them, where Cj only takes a value of zero or one. Thus 

\,k = ^i,k{^l, • • • ) Cj) 
= X! Mn{ci,--- ,Ci,Ci+i,--- ,Ck)Pp{uj{bi+i) = Ci+i,--- ,u{bk) = Ck) 

Ci+l,---,Ck 

Mn{ci,--- ,Ci-i,4,Ci+i,--- ,Ck)Pp{uj{bi) = c[,u}{bi+i) = Ci+1,--- ,u}{bk) = c^) 
= X] [M„(ci, • • • ,Ci, • • • ,Cfe) - M„(ci, • • • ,Ci_i,c-,Ci+i, • • • ,Cfe)] 

■Pp{u{bi) = c[,u{bi+i) = Ci+i,--- ,u{bk) = Cfc), (2.11) 
where the sum takes over all possible values of Cj and c^. On Qn{bi), 

Mnici+i, • • • , Ci, • • • , Cfc) - M„(cj+i, ■ ■ ■ , c- = 0, • • • , Cfc) = for Ci = or q = 1, 

Mn{Ci+i, • • • , Ci, • • • , Cfc) - MniCi+i, • • • , C- = 1, • • • , Cfc) = for C^ = Or Cj = 1. 

Thus, by (2.11), 

(l - Ig^ib.)) = 0, (2.12) 
where is the indicator of A- With this observation, 

= ^i,klg„{bi)- (2-13) 
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Note that on Gnih), 

[M„(ci, ■ ■ ■ ,Ci = 1, - ■ ■ ,Ck) - Mn{ci, ■ ■ • ,Ci_i,C- = 0,Ci+i, • • • ,Cfc)] = -1. 

so 

Ai,fe(ci, • • • ,Cj_i, l)/g„(6,) 

= ^6n(bi) [Af„(ci,- • • ,Ci = 1,- • • ,Cfc) - M„(ci,- • • ,Q_i,c- = 0,Ci+i,- • • ,Cfc)] 

•Pp(cj(6i) = 4 = 0,a;(6i+i) = Cj+i, • • ■ ,u{bk) = c^) 

+^Qn(bi) [Mnici,--- ,Ci = l,---,Cfe)-M„(ci,---,Ci_i,C- = l,Ci+i,---,Cfe)] 

•Pp(a;(6i) = c'i = l,w(6i+i) = q+i, • • ■ ,u}{bk) = c^) 

= ^GAb.) [^n(ci, • • • ,Ci = 1, • • • ,Cfe) - M„(C1, • • • ,Ci_l,c'j = 0,Ci+l, • • • ,Cfc)] 

■Pp{uj{bi) = c • = 0)Pp((:j(6i+i) = Q+i, • • • ,uj{bk) = Ck) 

= -{1 - p)Ig^^hi) 

Therefore, 

Ai,k{ci,---,Ci-i,l)Ig^(b,) = {l-pflg„ibiy (2-14) 
Similarly, note that on Qnih), 

[M„(ci, ■ ■ ■ ,Ci = 0, - ■ ■ ,Ck) - Mn{ci, ■ ■ • ,Ci_i,c'j = l,Ci+i, • • ■ ,Ck)] = 1. 

so 

Aj,fc(ci,---,Ci_i,0)/g„(fe,) 
= ^en(b») Y [^^n(ci,- • • ,Ci = 0, • • • ,Cfe) - M„(ci,- • • ,Ci_i,C- = 0,Ci+i,- • • ,Cfc)] 

•Pp(a;(6j) = c ■ = 0,uj{bi+i) = Cj+i, • • ■ ,u{bk) = Ck) 

+k„(b,) Y i^n{ci,-- • ,Ci = 0, • • • ,Cfe) - M„(ci,- • • ,Ci_i,C- = l,Ci+i,- • • ,Cfe)] 

•Pp(a;(5i) = c'i = l,uj{bi+i) = q+i, • • ■ ,u}{bk) = Cfc) 

= ^G^ib^) Y [^n(ci,---,Ci = 0,---,Cfe) -M„(ci,---,Ci_i,C- = l,Ci+i, •••,€/;)] 

■Pp{uj{bi) = 4 = l)Pp{uj{bi+i) = Cj+i, • • • ,uj{bk) = Ck) 

Therefore, 

Ai,fc(ci, • • • , Q-i, O)i'0„(5.) = P^Ig^ip^y (2.15) 
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Together with (2.13), (2.14) and (2.15), 

= ^p^f,fc(ci, • • ■ ,Ci-iA)I{ujib,)=i}Ign{b,) + Ep/\\j,{ci, ■ ■ • ,Ci_i,0)/{t^((,.)=o}^e„(6i) 
= Ep{l - pfl{oj{bi)=i}Ign{bi) + EpP^h^{bi)=o}Ign(bi)- 

Note that {a;(6j)} and Ig^{bi) independent, so 

EAlk = [(1 + -P)]^P^0„(6.) = [{l-p?P + pH^-p)]Pp{Gn{bi)). (2.16) 
By (2.10) and (2.16), 

k 

al{Mn) = Y,E,{Al,) = [{l-pfp + p\l-p)] PpiQnib)). (2.17) 

!=1 6G-Be{n) 

Therefore, if we divide the both sides of (2.17) by (2ra + 1)'^ and let n go to infinity, by (2.9), 
Jim a2(M„)/(2n + l)<^ = Ji^ [(1 - p)^^ + _ p)] ^ Pp(g„(6))/(2n + 1)<^ 

= -[{i-pfp+pH^-p)Wip)- 

So the theorem follows. 
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